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Abstract — In this paper, we show that there exists an arbitrary 
number of power allocation schemes that achieve capacity in 
systems operating in parallel channels comprised of single-input 
multiple-output (SIMO) Nakagami-ra fading subchannels when 
the number of degrees of freedom L (e.g., the number of 
receive antennas) tends to infinity. Statistical waterfilling - i.e., 
waterfilling using channel statistics rather than instantaneous 
channel knowledge - is one such scheme. We further prove that 
the convergence of statistical waterfilling to the optimal power 
loading scheme is at least 0(1/ (LlogL)), whereas convergence 
of other schemes is at worst 0(l/logL). To validate and 
demonstrate the practical use of our findings, we evaluate the 
mutual information of example SIMO parallel channels using 
simulations as well as new measured ultrawideband channel data. 

Index Terms — Capacity, power loading, diversity, statistical 
channel state information, Nakagami-m, UWB. 



I. Introduction 

Parallel channels are frequently encountered in mod- 
ern communication systems. Examples of such systems in- 
clude those using orthogonal frequency division multiplexed 
(OFDM) transmissions, time-division multiplexed (TDM) 
transmissions, frequency-hop spread-spectrum (FH-SS), and 
multiple-input multiple-output (MIMO) systems that employ 
eigenmode transmission techniques. 

In recent years, the design of subchannel power allocation 
strategies has become a popular research topic for systems 
operating in parallel channels. Power loading schemes based 
on having perfect channel state information at the transmitter 
(CSIT) have been developed to maximize spectral efficiency 
(T) and system throughput J2, as well as to minimize bit- 
error rate 0, flU and transmit power [5|. Solutions that rely 
on partial CSIT have also been devised (see, e.g., J6)). 

Here, we investigate power loading in parallel channels 
using only statistical CSIT (SCSIT). Several researchers have 
reported notable results on the use of SCSIT in the literature. 
For example, power loading based on SCSIT to minimize 
bit-error rate was studied in Q, while in (8), the authors 
developed power and bit allocation algorithms to maximize the 
spectral efficiency and minimize the power consumption. More 
recently, in |9|, the authors investigated transmission in jointly 
correlated MIMO channels, whereby the transmit vector is 
conveyed over the eigenmodes of the transmit correlation 
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matrix, and power allocation is performed to maximize an 
upper bound on the ergodic capacity. It was found that, 
for parallel channels in particular, statistical waterfilling - 
i.e., waterfilling using the mean channel gains instead of the 
instantaneous gains - optimized this bound [9, cf. (67)]. 

In this paper, we study power loading in parallel chan- 
nels comprised of single-input multiple-output (SIMO) sub- 
channels where the number of degrees of freedom of each 
subchannel (e.g., the number of receive antennas) is denoted 
by L. We maintain generality in our analysis by assuming 
the single-input single-output (SISO) channels that consti- 
tute each SIMO subchannel adhere to a Nakagami-m fading 
model. In contrast to J9), we investigate conditions under 
which we can determine the power loading schemes that 
maximize the exact ergodic capacity of such a parallel channel 
rather than an upper bound. To this end, we show that as 
L — > oo, there are arbitrarily many such schemes, where 
statistical waterfilling is one such approach. Furthermore, we 
demonstrate that convergence to the optimal power allocation 
scheme is 0(1/ (LlogL)) for statistical waterfilling, whereas 
the convergence of other asymptotically optimal schemes is 
at worst 0(l/logL). Finally, we validate and demonstrate 
the practical use of our findings through a study of SIMO 
parallel channels using simulations as well as new measured 
ultrawideband (UWB) channel data. 

The rest of the paper is organized as follows. Bounds on the 
capacity of a parallel channel are given in Section [ID These 
bounds are used to develop the main results on the asymptotic 
optimality and convergence properties of statistical waterfilling 
and other power loading strategies in SIMO systems in Section 
Hill Measured and simulated channel data are used to measure 
performance in Section HVl and conclusions are then drawn in 
Section |V] 

II. Bounds on Capacity 

Consider a parallel channel comprised of N SIMO sub- 
channels, each with L degrees of freedom^ The mutual 
informatioijl of the nth subchannel, assuming full channel 
state information is available at the receiver, is given by [ 1 1 

I ln (P n ) = log + ^ 7 „ 

where j n = ||h„|| 2 with h„ e C L being a vector of complex 
channel coefficients, P n is the power transmitted on this 

The restriction that all subchannels have L degrees of freedom is made 
for simplicity of presentation; in fact, the main result of this paper holds even 
if L varies across the subchannels. 

2 Mutual information and capacity are defined in units of nats. 
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subchannel, and Nq is the variance of the zero-mean, additive 
white Gaussian noise on that subchannel It follows that the 
ergodic capacity of the parallel channel is given by 



C 



sup 

Pn>0, 



■ N 
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„ n— 1 



E [I 7b (Pn 



(1) 



Since only SCSIT is available, the powers {P„} cannot be 
functions of {7„}, but instead rely upon the statistics of these 
channel gains. We define the set {P*} as the set of optimal 
powers, i.e., the set that yields the supremum given above. 

To develop our main results, we will require an upper and a 
lower bound on C. An upper bound is easily calculated using 
Jensen's inequality. Let P* denote the power loading strategy 
based on statistical waterfilling, which is given by 

where (x) + = max{0,x}, fj, n = E [7J is the mean of the 
nth subchannel gain, and v is chosen to satisfy P* = P. 
Now, we can write 



N 



JV 
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(Pn) 



(3) 



where the first inequality follows from the concavity of / and 
the second results from the standard waterfilling solution. 

Clearly, {P£} is not the only set that satisfies this bound; 
it is simply the set that maximizes the bound. It will be 
convenient to define another set of powers {P,' } that satisfies 

N N 

C<£WO<£W^) (4) 

n— 1 n— 1 

such that P n = P- Note, we do not explicity define {P^}. 
For the lower bound, we begin with the trivial inequality 



c> > p[/ 7 „ (i*; 



n=l 



(5) 



Applying Markov's inequality yields the more useful bound 



n=l " 



r 



(6) 



where JF n (x) = Pr (7„ > x) is the complementary cumula- 
tive distribution function (CCDF) of the nth subchannel gain 
and a n > 0. Note that a sufficient (but not necessary) condition 
for this bound to hold is P,' = P* for all n. 

III. Power Loading in Diversity Channels 

We now present our main results. We consider the case 
where each subchannel experiences Nakagami-m fading. The 
constituent single-input single-output (SISO) channels corre- 
sponding to the nth SIMO subchannel are statistically inde- 
pendent, and the £th SISO channel gain, denoted by |ft.„,f| 2 , 
is gamma distributed with scale parameter 8 n and shape 

3 We assume without loss of generality that the power spectral densities of 
the noise processes on all subchannels are identical. 



parameter m n . Note that we assume the scale and shape 
parameters are identical for all SISO channels that constitute 
the nth subchannel. Thus, the channel gain 7„ = ||h„|| 2 is a 
gamma distributed variate with scale parameter 8 n and shape 
parameter m n L. The density function of j n is given by 



fnh) 



-,ro„L-lg-5ir 



, 7> 0. 



(7) 



d^ L T(m n Ly 

The mean of the nth subchannel gain in this case is given by 



Mn 



,m n L 



(8) 



Now we can study the capacity as the number of degrees 
of freedom L grows large, which leads to the following 
proposition. 

Proposition 1: Consider a parallel channel with subchannel 
gains distributed according to (0. In the limit of large L, there 
are arbitrarily many asymptotically optimal power loading 
strategies, one of which is the statistical waterfilling solution 
defined by ©. 

Proof: We prove the proposition by showing that the 
upper and lower bounds given by 



N 



N 



(9) 



n=l 



n=l 



are asymptotically equivalent (as L —> 00) for any set {P,'} 
that satisfies (O for L > Lq. Note that we do not explicitly 
define {P n }, although it is clear that one possible definition 
is P' n = P,j for all n since this definition satisfies (O for 
all L. We first consider the upper and lower bounds on the 
capacity of the nth subchannel. We assume that a nonzero 
power is allocated for transmission on this subchannel; this is 
a valid assumption since if the converse were true, the upper 
and lower bounds on capacity would both be zero. The CCDF 
of the nth subchannel gain is 



(10) 



F n (x)=Q \ m n L, — 



where Q (a, x) = T (a, x) /Y (a) is the normalized upper 
incomplete gamma function. Now, letting /3 n = P^9 n m n /No 
where P' n > 0, we can write the ratio of the lower and upper 
bounds given by (|6]l and ©, respectively, as a function of L: 



tl = 



a n Q (m n L, fi n 1 m n (e a " - 1)) 
log (1 + p n L) 



(ID 



For every LeN and n <E {1,2,..., A^}, we can calculate j3 n . 
Thus, we can choose to define 



a n = log (1 + a/3 n L) > 

where < a < 1. It follows that 
\og{l + a/3 n L) 



log 



(m n L, am n L) . 



(12) 



(13) 



Using HQ 5.11.3] and flj] 8.11.6], and noting that < 
< 1 for < a < 1, we can expand at 00 to 



a c 



l-a 
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logarithmic term 

x m n L 
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gamma function term 



Thus, tl — > 1 as L — > oo. Since this relation holds for any 
n, the upper and lower bounds given in Q are asymptotically 
equivalent, and {P„} is an optimal set of powers. Finally, it 
follows from continuity that any set {P r ' + e„} with e n € R, 
which satisfies (O, is also optimal. □ 

This result can be understood intuitively by noting that the 
mean channel gains increase monotonically with the number 
of degrees of freedom L. Thus, the effective SNR on a given 
subchannel grows without bound as L — > oo, suggesting that 
an equal power allocation strategy is optimal. Many power 
loading strategies satisfy this condition, including traditional 
waterfilling and statistical waterfilling. In fact, this result is 
also related to the notion of "channel hardening", which has 
been studied in the information theoretic literature. In par- 
ticular, our result is corroborated by related findings detailed 
recently in JHJ cf. §IV]. 

Proposition [T] is stronger than it may first appear since it 
implies that one may theoretically determine the number of 
diversity branches such that, when an asymptotically optimal 
power loading strategy is employed, the ergodic capacity can 
be approached to an arbitrarily close degree. More exactly, for 
arbitrarily small e > 0, there exists some positive integer Lq 
such that when L > Lq we have \C — C'\ < e, where C 
denotes the rate obtained by employing the suboptimal power 
allocation scheme {P^}. 

Although we have demonstrated the asymptotic optimality 
of an arbitrary number of power loading schemes, including 
statistical waterfilling, it is essential from a practical viewpoint 
that we understand the convergence behavior of this result. 
Convergence can be studied by observing the rate at which the 
upper and lower bounds on capacity, given by ||9), approach 
the limit in L. This leads to our second main result. 

Proposition 2: Statistical waterfilling converges to the opti- 
mal power loading strategy like 0(1/ (LlogL)), whereas for 
more general asymptotically optimal power loading strategies, 
convergence is, at worst, O (1/logL). 

Proof: The proof of the first part follows by substituting 
P' n = P£ in the expression for Tl given by ( TT31 and expanding 
at L = oo. The second part follows from the proof of 
Proposition [TJ. □ 

From this proposition, we see that although there are arbi- 
trarily many asymptotically optimal power loading strategies, 
convergence to optimality is guaranteed to be significantly 
faster for statistical waterfilling than for general schemes. 
Consequently, statistical waterfilling is not only easily imple- 
mented and intuitive, but is a practical solution for systems 
with moderate levels of inherent diversity. 
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Fig. 1. Interconnection of sounding equipment. 



IV. Numerical Evaluation 

One example of a parallel channel with subchannels that 
experience unequal mean fading gains can be found in OFDM- 
based UWB systems. In this section, we evaluate our findings 
by considering both measured and simulated UWB channels. 

A. Measured Channels 

A channel measurement campaign was conducted using a 
state-of-the-art 2x4 time-domain multi-antenna UWB channel 
sounder. The sounder, manufactured by MEDAV, interrogates 
the propagation channel by using trains of pseudo noise (PN) 
sequences of 4095 pulses or chips [13]. These are generated in 
baseband at a clock rate of 6.95 GHz and later up-converted, 
using the same clock, to cover the bandwidth from approxi- 
mately 3.5 to 10.5 GHz. In turn, the receiver down-converts 
captured signals, does a periodic sub-sampling of them, and 
uses a phase shifter to allow the sampling of complex channel 
impulse responses (CIRs) in the time domain. 

Fig. Q] shows the configuration of the sounding equipment 
for measurements. Test signals are generated and bandpass 
filtered to avoid out-of-band emissions. The test signals are 
then amplified, transferred to a biconical antenna for radiation 
and, after travelling through the propagation environment, 
received by a similar antenna connected to a bandpass filter 
lfl4l . fTBI . After correcting the gain of the incoming signals 
with an automatic gain control (AGC) unit, the receiver uses an 
analogue-to-digital converter (ADC) and performs a matched 
filtering of the data, with the known PN sequence, in a digital 
signal processing (DSP) unit. Before sounding and after a 
warming-up period, the system response, phase imbalance and 
crosstalk are characterized using cabled or open connections. 
These measured parameters are then used for calibration to 
leave only the combined response of the antennas and the 
environment on the data recordings. 

UWB (2 x 4) MIMO channel measurements were conducted 
in an open-plan modern office in the central area of Bristol, 
UK. The sounding environment had typical scattering objects 
such as personal computers, liquid crystal displays (LCDs), 
non-metallic cubicle partitions, desks and metallic cabinets. 
The ceiling of the office was made of perforated metallic tiles 
and the floor supported by a metallic structure covered with 
non-metallic material. The dimensions of the sounded office 
are: 12.74 m wide, 30.84 m long and 2.39 m high. 

For this work, each recorded CIR results from averaging 
256 captured CIRs in hardware. Under these conditions, the 
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Fig. 2. Illustration of sounded environment and antenna locations. 



observation time of the system for a recorded CIR is 155 ms. 
The antennas were mounted on fibreglass masts at 1.3 m 
from the floor. The transmit antenna mast was attached to the 
movable part of an x — y automated positioning system 11611 . 
The positioners and the sounder were remotely controlled to 
prevent human intervention in the area of measurements. In 
parallel, a spectrum analyzer, connected to a biconical antenna 
and a low noise amplifier, periodically scanned the spectrum of 
interest. This information was was used to check the "health" 
of the test signal and the presence of interfering signals. 

Fig. |2] shows a diagram of the sounded environment. Two 
specific locations were selected to capture line-of-sight (LOS) 
and non-line-of-sight (NLOS) data. These locations, T1 (LOS) 
and T2 (NLOS), are grids of x—y points in which two transmit 
antennas, 25 cm apart, were displaced at distances of 3 cm. 
The receive antenna mast was positioned at a fixed location 
(Rx), and a four-element linear antenna array was formed by 
securing the first antenna and the following ones at 3, 6 and 
12 cm from it. Each measurement grid had 441 (21 x 21) 
points, so a total of 3528 (441 x 8) CIRs were recorded for 
each location T. In general, the distance between transmit and 
receive antennas ranges from approximately 7.02 m to 6.51 m 
for T1 and from 6.41 m to 5.77 m for 12. Note that NLOS 
conditions were achieved by locating the transmit antenna in 
such a way that a concrete column was always shadowing it 
from the receive one. 

SIMO channels are formed by extracting the CIRs measured 
from the first transmit antenna to the receive antennas spaced 
apart by 12 cm distance. Each of these CIR is converted into 
the frequency domain, by using a discrete Fourier transform 
(DFT), and we retain the channel frequency response charac- 
teristics for the 5-to-6 GHz band. With a frequency spacing 
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Fig. 3. Plot of the power of three consecutive channel frequency response 
measurements in the 5.3-IO-5.5 GHz band. 
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Fig. 4. Mean channel gains for LOS and NLOS measured channels in the 
5-to-6 GHz band. 



of roughly 1.7 MHz, this amounts to 588 frequency samples, 
i.e., channel frequency response coefficients. A power plot of 
three consecutive NLOS measurements in the 5.3-to-5.5 GHz 
band is illustrated in Fig. [3] From this figure, we can see that 
the measurements are fairly independent of one another. 

The channel measurements were normalized such that the 
average mean channel gain (in frequency) is one. Fig.|4]depicts 
the mean channel gain (averaged over the 441 available snap- 
shots) for the LOS and NLOS channels. It is certainly clear 
from this figure that the mean fading gains vary considerably 
with frequency. Thus, one would expect an unequal power 
loading strategy based on the statistics of these channels to 
perform better than a balanced power allocation scheme. 

In Fig. [5] bounds on the capacity, defined by (0) and 
0, are plotted for the measured LOS and NLOS channels. 
Additionally, the mutual information for a channel where a 
balanced power allocation is used is illustrated. The lower 
bound and the balanced power rate are both averaged over the 
441 snapshots of measured data to emulate the expectation 
in the rate expression. Since a large range of SNR values 
are considered in these graphs, it is beneficial, for ease of 
comparison, to normalize these results with respect to the 
capacity of a parallel additive white Gaussian noise (AWGN) 
channel. It is evident from Fig. [5] that statistical waterfilling 
is capable of providing significant gains at low SNR, which 
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Fig. 5. Normalized mutual information (w.r.t. a parallel AWGN channel) 
and bounds for measured LOS channels (a) and NLOS channels (b). 
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Fig. 6. MPE for measured LOS and NLOS channels. 



exemplifies typical operating conditions in wideband systems 
with a strict power budget, such as UWB. 

It is also useful to examine the maximum percent error 
(MPE) of the capacity bounds, which is defined as MPE = 
100% x (Cub — Cub) /Clb where Cub and Cub are the 
upper and lower bounds given by (0 and (0, respectively. This 
metric quantifies the deviation of the statistical waterfilling 
power allocation strategy from the optimal power allocation 
since the power-optimal rate lies between Cub and Cub- The 
MPE is illustrated in Fig. [6] which shows that even for the 
practical situation of L = 4, the difference in the two bounds 
is relatively small, from which we deduce that statistical 
waterfilling is a good pragmatic power loading strategy. 

B. Simulated Channels 

In the measured scenario, we are limited to studying the 
mutual information of a SIMO channel with L < 4 receive 
antennas. Thus, we employ simulations to observe the perfor- 
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Fig. 7. MPE for measured LOS and NLOS channels. 



mance of various power loading strategies as a function of 
L. To this end, we consider a system operating in the 5-to-6 
GHz band as above, but where the mean channel gain decays 
with frequency like / 3 , which is typical in rich scattering 
environments. 

We are particularly interested in the MPE of the upper and 
lower capacity bounds, and plot this metric in Fig. [7] as a 
function of L for SNR values of —10 dB and 5 dB. This figure 
confirms that statistical waterfilling converges to the optimal 
power loading strategy quickly, achieving an MPE of less than 
5% for L > 4, thus making this technique suitable for many 
practical systems. 

V. Conclusions 

In this paper, we showed that there exist an arbitrary number 
of power loading strategies that achieve capacity for parallel 
channels comprised of SIMO subchannels in the limit of a 
large number of degrees of freedom L. Statistical waterfilling 
was shown to be one such strategy, which in fact is guaran- 
teed to converge to the optimal solution much quicker than 
other asymptotically optimal solutions. We demonstrated the 
practicality of our results through a study of measured and 
simulated UWB channels. 
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